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In two recent papers, arXiv:0909.1559 and arXiv:091 1.4122, it has been demonstrated that one can ob- 
tain quantum corrected Boltzmann kinetic equations for leptogenesis using a top-down approach based on the 
Schwinger-Keldysh/Kadanoff-Baym formalism. These "Boltzmann-like" equations are similar to the ones ob- 
tained in the conventional bottom-up approach but differ in important details. In particular there is a discrepancy 
between the CP-violating parameter obtained in the first-principle derivation and in the framework of thermal 
field theory. Here we demonstrate that the two approaches can be reconciled if causal n-point functions are used 
in the thermal field theory approach. The new result for the medium correction to the CP- violating parameter is 
qualitatively different from the conventional one. The analogy to a toy model considered earlier enables us to 
write down consistent quantum corrected Boltzmann equations for thermal leptogenesis in the SM+3vr which 
include quantum statistical terms and medium corrected expressions for the CP-violating parameter. 
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I. INTRODUCTION 
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To calculate the baryon asymmetry generated during the epoch of leptogenesis [1] in the standard model 
extended by three right-handed neutrinos (SM+3^r) and its extensions one usually uses standard Boltzmann 
kinetic equations. The collision terms (and in particular the CP-violating parameters) in this equations are 
computed in vacuum in the in-out formalism [2, 3] and do not take into account effects induced by the hot 
medium of the early universe. Such effects can be consistently taken into account in a top-down approach 
based on the Schwinger-Keldysh/Kadanoff-Baym formalism. In [4, 5] we have applied it to a simple toy 
model of leptogenesis and derived a new (quantum corrected) form of the Boltzmann equations, which 
includes quantum statistical factors and takes the medium effects into account. We have found that the 
medium corrections to the CP-violating parameter e depend only linearly on the one-particle distribution 
. functions (see also [6-8]). In the analysis based on finite temperature field theory for the phenomenological 

scenario of thermal leptogenesis [2, 3, 9] and for GUT baryogenesis [10] the medium corrections to the 
. CP- violating parameter e depend quadratically on the distribution functions. 

This discrepancy has been noted in the context of leptogenesis in [4] for the vertex contribution to the 
CP-violating parameter and later in [5] for the self-energy contribution. Here, we use a finite temperature 
equivalent of the Cutkosky cutting rules [11-14] to derive thermal corrections to the expression for the 
imaginary part of the three-point vertex function and the self-energy loop and to calculate the corresponding 
medium-corrected CP-violating parameters. We show that the discrepancy is due to an ambiguity in the 
real-time (RTF) formulation of thermal quantum field theory and disappears if one considers retarded or 
advanced n-point functions. In the framework of the toy model this has been demonstrated recently in [15]. 
Together with the new form of the Boltzmann equation derived in [4, 5] this puts us in the position to write 
down quantum corrected Boltzmann equations for the phenomenological scenario of thermal leptogenesis 
which consistently include the medium corrected CP-violating parameter and quantum statistical terms. 

In section II, we introduce our notations for the CP-violating parameters and the thermal field theory for- 
malism. Then, in section III we review the conventional calculation of the thermal corrections, and in 
section IV we demonstrate how to reconcile them with the recent results from nonequilibrium field theory. 
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Finally, in section V we present the quantum corrected Boltzmann equations taking medium corrections 
into account. 



II. CP-VIOLATING PARAMETER AND THERMAL FIELD THEORY 



In the phenomenological scenario of thermal leptogenesis as well as in the toy model, the matter-antimatter 
asymmetry is generated by the decay of a heavy species. In both cases the CP violation in this decay is 
caused by interference between the tree level and the one-loop diagrams, see fig. 1. In the phenomenological 




FIG. 1: Tree level and one-loop contributions to the heavy Majorana neutrino decay ipi —t ctf3. The asymmetry, at 
lowest order, is due to the interference of these contributions. 

scenario ipi — N{ are heavy Majorana neutrinos which decay via Yukawa interactions C = h a i£ a Ni <p+h.c. 
into leptons a — I and Higgs @ = <p or their anti-particles. In the toy model ipi is a heavy real scalar particle 
which decays via Yukawa interactions C = — ^ipibb+h.c. into two light scalars a = f3 = b or the conjugate 
b. The CP-violating parameter for the decay of ipi is defined as 

i i ' i J-* , -p _ ) \ l J 

where rjy 4 _^ includes a sum over flavour indices and loop-internal Majorana neutrino generations in the 
case of the phenomenological scenario: r^-^ = J2 a j ^Ni^e a <t> (we do not consider flavor effects here). 

If the tree level and one-loop contributions are written as Ao-4o and Ai^i, respectively, where all coupling 
constants are absorbed in Aom, the CP-violating parameter becomes at lowest order: 

\\ Ao + XiAil 2 ~ \X*A + XIA^ 2 ^ JmjAgA^Iml^A} 
e ' lAo^lo + AiAl' + IA^o + AI^I 2 ~ |A | 2 |A| 2 

where the sum over lepton flavour and Majorana neutrino generation indices is again implicit. In the case 
of thermal leptogenesis this leads to 

where q and k denote the four-momenta of the Majorana neutrino and the lepton respectively, see fig. 3, and 
for the toy model to 



e i = -2\g j \ 2 Im[ ^ )lm{^i} , Mi, » = 1,2. (4) 
\9i 9j J 

This means that one needs to compute the imaginary (absorptive) part of the vertex and the self-energy loop 
contributions Im {„4q„4i}. In vacuum this can be done conveniently with help of the Cutkosky cutting rules 
[16-18]. In thermal quantum field theory these can be generalized in order to take into account interactions 
of internal lines in the loops with the background medium [11, 13, 14]. In the real-time formalism of 
thermal quantum field theory two types of fields, termed type-1 and type-2 fields, are introduced in order to 
avoid pathological singularities [19]. Vertices can be of either type, differing only by a relative minus sign. 



9i9i 
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We denote them by g 1 = —ig and g 2 = +ig for a generic coupling 1 g. The propagators connecting the 
different types of vertices can be considered as components of a 2 x 2 propagator matrix 2 

For a scalar particle b the components are 

A b (p) = D b (p), A±(p)=Z?±(p). (5) 
For a fermion / the components are 

A f {p) = ^-p + m f )D f (p), Af(p) = (7-p + m f )Df(p). (6) 
For brevity we have defined 

D a (p) = - l - T — - 27rU a ' eq (p)S(p 2 - ml) , 

Dt(p) = 2ir [Q(±p ) - U a ' eq (p)} S( P 2 - m 2 a ) , (7) 

where £ a = +1 for fermions and £ a = —1 for bosons. Here, we denote by f b ' eq (p) and f^ ,eq (p) the 
equilibrium distribution function for bosons and fermions, respectively, given by 

f a ' eq (p) = [exp(/3\p^\)+^}- 1 . (8) 

They are functions of the Lorentz invariant product p^U 11 of the particles' four-momentum and the four- 
velocity U of the plasma in a general frame. In the rest-frame of the plasma, U = (1, 0, 0, 0), we obtain the 
standard form which depends on po- In tne following we assume that it is sufficient to replace the different 
propagators in our toy model and the phenomenological theory by their thermal field theory equivalents 
given in eqn. (5). This approach has been followed in previous works for the baryogenesis and leptogenesis 
scenarios [2, 3, 9, 10]. We ignore further thermal effects, such as thermal corrections to the masses and 
wave function renormalization here. 

Denoting vertices attached to external lines by Xi and those attached to internal lines only by Zj we can 
formally denote an amputated n-point graph by F(xi, . . . ,x n ; Zj). Here we assume that F is given in 
momentum space, writing the position space coordinates in order to identify the individual vertices. The 
contribution of this graph to the amplitude is —iF(xi, . . . , x n ; Zj). 

Physical amplitudes involve a sum over possible combinations of types of internal vertices: 

J~ {x\ , . . . , x n , Zj) ^ F(x\ , . . . , x n , Zj) . 

type Zj 

For external vertices of fixed type it has been shown [11, 12] that this sum is equivalent to a sum over all 
possible "circlings" of the internal vertices: 3 

T(xi, . . . ,x n ;zj) = ^2 F^(xi,...,x n ;zj). (9) 

circling Zj 

F> and F< with "circled" vertices represent graphs computed using the set of rules, given in fig. 2. 
These differ for the computation of F> and F< by interchange of the A + and A~ propagators. In 



1 At the end of the calculation of Im {^IJ^li } we set g = 1, since the physical coupling constants have been factored out into Ao 
and Ai. 

2 In [2] and elsewhere resummed propagators have been used in this place to prevent the appearance of singularities. Since we are 
mainly interested in the structure of the thermal corrections we stick to the free thermal propagators here. 

3 The historic origin of this formula was that the external fields where considered to be all of type 1 (physical). 



4 



F^(xi, . . . , x n ; Zj) we explicitly denote circling of a vertex a as F^(xi, . . . , x Q , . . . , x n ; Zj). Note that 
the two ways of defining T in terms of F > and _F< in eqn. (9) are in agreement only if the Kubo-Martin- 
Schwinger (KMS) boundary condition, 

A-(p) = -Z a e-^ u AZ(p), (10) 

is satisfied. This is the case in thermal equilibrium. 




= A+(p) 



G =A-(p), G O = A*(p) 

P P 

FIG. 2: Circling rules for a generic theory used for the computation of F> in momentum space. The rules for the 
computation of F < differ by interchange of A + (p) and A~ (p). The A propagators connecting circled and uncircled 
vertices may be interpreted as cut propagators. In vacuum they correspond to the cut propagators in the Cutkosky rules. 

From T we can then compute Iiil{„4q*4i} as 4 

Imfy^A} = -\m{^—\ . (11) 

1 515233 J 

where gx, <?2 and g% stand for the generic couplings associated with the three vertices in the one-loop 
diagrams at Fig. 2. 

III. PHYSICAL AND GHOST FIELDS 

In this section we briefly review the conventional calculation of the CP-violating parameter in real-time 
thermal field theory. However, we use a notation that is helpful to understand the ambiguities emerging 
there, and that can be more easily compared to the results from non-equilibrium field theory. 

An obvious problem with the real-time formulation for the computation of n-point functions is that there 
are in general 2" such functions which differ in the types of the external vertices. Historically the correct 
function was considered to be the one with all external vertices of type-1 (physical). In this case eqn. (9) 
leads to the following formula for the imaginary part of a graph's contribution to the amplitude: 

Im{£ _1 .F(l,...,l;z,-)} = - F-^(x 1 , . . . ,x n ; , (12) 

circling (x^), Zj 

where the sum includes all possible circlings of the internal vertices Zj but only those circlings of external 
vertices Xi which include both, circled and uncircled vertices (indicated by the brackets around Xi). 

The six diagrams contributing to the imaginary part of the three-point vertex function according to eqn. (12) 
are shown in fig. 4. The circlings contributing to the self-energy part are shown in fig. 5. 



4 In the phenomenological scenario the Feynman rules for Majorana neutrinos include spinors, charge conjugation and projection 
operators which we assume to be included in _F> . 
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f q+ I— k 




FIG. 3: Momentum flow in the vertex and the self-energy loop. 




(d) 



(e) 



(f) 



FIG. 4: Circlings contributing to Im {i~ 1 T(l, 1, 1)} for the vertex loop. At one-loop level the circlings can be inter- 
preted as cuts, as indicated, by the lines separating circled from uncircled regions [11]. The contributions from diagrams 
involving cuts through the X2 -X3 line are suppressed relative to the others in the hierarchical limit. 




(a) 



(b) 




(c) 




(d) 



FIG. 5: Circlings contributing to Im {i~ \F(1, 1; z)} for the self-energy loop. The graphs (b) and (c) vanish since ipi 
and tpj cannot be on-shell simultaneously. Note that we consider only the diagrams with ipi in the external and ipj in 
the internal line (i 7^ j) because these are the only ones which contribute to £j. 



The contributions which correspond to cuts through the ipj line are suppressed in the hierarchical limit. If 
they are neglected the application of this circling formula leads for the toy model to the result 



8tt M( 

for the vertex contribution and 

2 



919 3 



S,th 



\9^(9i9 3 



1 



4tt Mj/M? + |(1 + cos 0i) 



^{1 + /If + fk q + 2/If /Jf } (14) 



(13) 



Hi-"" V '/;//, / .\/f Mf .1 4n 
for the self-energy contribution. The distribution functions are to be evaluated for the energies E\ and E2 
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given by 

El2 = -[Ep T | q| (sin t cos tpi cos 6' + cos 0, sin 8')] , (15) 

where 9i and ipi are elements of the solid angle Hi and the angle 6' is given in the limit of massless decay 
products by sin<5' = (|p| — |k|)/|q|. The dots in eqn. (13) represent further terms in f^ j ' eq which are ne- 
glected. Equivalently these results can be derived directly using only the 11 components of the propagators, 
because the vertex and the self-energy loop do not include internal vertices. Very similar results are known 
for the phenomenological scenario which can be obtained using the propagators in eqn. (6) for fermions for 
the Majorana neutrinos and leptons in the loops and eqn. (5) for the Higgs bosons [2, 9]. For the dependence 
on the distribution functions one obtains then the quadratic form 

i-/if+/ir-2/if/ir. (i6) 

The results obtained from non-equilibrium field theory in [4, 5] differ from eqns. (13) and (14). The non- 
equilibrium results feature a different dependence on the distribution functions, 

i + fk+fk + 2 fhfh -> i + + d7) 

Note that the top-down results are valid even if f b is not an equilibrium distribution (f b ~ f b must hold, 
however) and that the dependence is linear in the distribution function in contrast to eqns. (13), (14). The 
latter property contradicts the result derived from thermal quantum field theory . 

In the phenomenological model, an analogous replacement leads to a particularly important discrepancy. 
Indeed, eqn. (16) would imply a cancellation of the leading effects since fp eq — f^ eq — 2fp eq fp ,eq . The 
remaining effect is, in this case, entirely due to the fact that different energies enter the distribution functions 
of leptons and Higgs particles in eqn. (16). Since E\ — E% ~ |q|, this effect vanishes when the velocity 
of the Majorana neutrino in the medium rest-frame, {ql/E^ 1 , becomes small. Therefore it is important to 
check wether a replacement of the form of eqn. (17) does also occur in the phenomenological scenario. This 
will be investigated in the next section. 



IV. CAUSAL N-POINT FUNCTIONS 



We will now see how the finite temperature field theory approach can be reconciled with the results derived 
from non-equilibrium quantum field theory. In [21-23] it was shown that the combination 

3~ r/a(^i ' ■ * ■ ' z j ) ^ ^ E^z [x\] . . . x a , . . . , x n ] Zj) , (18) 

circling Xi : Zj 

referred to as the retarded (advanced) product, has the distinguishing property that the time component 
(x a )o is singled out as being the largest (smallest). This becomes clear when we consider the so-called 
largest (smallest) time equation 

F^(xi,.. . ,x a> ... ,x„) + F^(xi,.. . .,£„) = 0, if (x a ) largest/smallest, (19) 

which implies pairwise cancellation of the terms in eqn. (18) if any external vertex Xi with i ^ a has the 
largest (smallest) time component. It has been realized that such causal products appear in Boltzmann 
equations in different cases, see for example [22, 24]. Furthermore, it has been shown that the causal 



1 Note that if in eqn. (14) the term quadratic in the distribution functions was absent, then by redefining the integration variable ipi 
we could write the energies E± and E2 in the form E\ t 2 = \ {Eq 1 + |q|(sin 0; cos ipi cos 8' =p cos 81 sin <5')] , which was used 
in [20]. 
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products agree with the results of the calculation in imaginary-time formalism analytically continued to 
real energies, at least in a few examples including the self-energy loop and the three-point vertex. 

The imaginary part of the causal sum was shown in [22] to obey 



Im{i l T%) A {x 1 ,.. 



T 2 



not all 

E E 

circling x% circling Zj 



, X n J Zj ) } — 

Im|r 1 F > (xi,... 



i X n \ Zj) 



F. 



<\X\, ■ ■ ■ , x Q 



(20) 



where "not all" means that not all Xi should be circled at the same time and the imaginary part is taken of 
the causal product in momentum space. Here, the vertex x a with largest or smallest time is always circled. 

We can now compute the imaginary part of the advanced product Im{z _1 J r ^ 1 ^(xi, x%, £3)} for the three- 
point vertex with smallest time component {xi)q of the decaying particle. The relevant circlings are shown 
in fig. 6. As before the contributions fig. 6(b) and (c) are suppressed due to the cut through the tpj propagator 
line. 

X2 




3* 



(a) 



(b) 



(c) 



FIG. 6: Circlings contributing to Im{i 1 J r ^ (xi,X2, £3)} for the vertex loop. The advanced three-point function in- 
volves a difference of F > and F< contributions which differ by the replacement A <-» A T in the circling rules. Since 
the finite temperature contributions (terms proportional to f eq ) to the latter are the same, all contributions quadratic in 
the distribution functions cancel. 



We compute the remaining contribution from fig. 6(a): 

Im{^ 1 J r ^ 1) (a;i,X2,a;3)} =-Im{i~ l F > (x_ 1 , x 2 , x 3 ) - i~ 1 F < (x 1 , x 2 , x 3 )} 



dn r 



(27T) 



( + i gi )(-ig 2 )(-ig 3 )Dp(q + l)D^(q + I - k)D+(l)- 



- (+m)(-m)(~m)D^(q + i)D^( q + 1 - k)D-(i)] s] 



(21) 



where we take F^ to include the spinors and charge conjugation operators C as well as projection operators 
Pr, Pl associated with the vertices (for the Majorana neutrino interactions). This leads to the trace part 
denoted by S [2, 25, 26]. In the massless lepton and Higgs limit: 

S=Y, [ui{k)P L u N MY[Mk)PL(j -(q + l-k) + M^C-'Pl^ ■ l)P R Cu Ni {qj\ 

spins 



=Tr[(7-g-M i )(7.fc)Mj( 7 .0Pfi 



2MiMjk-l 



(22) 



and S = 1 for the toy model. 

It turns out that the pole of the i/jj propagator does not lie in the loop integration region, so we can drop the 
it prescription. We then get (the upper and lower signs correspond to the toy model and phenomenological 
scenario respectively) 
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lm{i 1 J r< x\xi,x 2 ,x 3 )} 



-— Im 
2 



dH 



S((q + l) 2 -ml)6(l 2 -m 2 a )- 



(2tt) 2 u * ; pr y °" (q + l-k) 2 - M. 

{ [e(-(q + l ))Q(l ) - e(q + lo)Q(-lo) 



± 



e(-( qo + i ))-e(q + i ) fr q + 



(23) 



which becomes 



-if W? S{{ " + "' - ^ ~ ^ (, + < - 1)' - {' * ^ + ! "») S - <24> 

Performing the integration over dlo, this indeed leads to the result for the CP-violating parameter 
obtained in the top-down approach with correct dependence on the distribution functions eqn. (17). For the 
phenomenological scenario we obtain in the limit of massless lepton and Higgs: 



v.th = J_ Im {( h * h )i 3 } M 3 



16tt 



where E\p are given by eqn. (15). In the zero temperature limit this reduces to the well-known result 

8tt 



Mi 



1 — COS 6i 



4tt Mf/M? + + cos 0i) 



(25) 



Va 



-E 



3& 



with 



1 - (1 + x)\n 



1 + x 



(26) 



The same computation can be performed for the self-energy loop. The possible circlings are shown in fig. 7: 





(a) 



(b) 



FIG. 7: Circlings contributing to Im{i 1 J r ^(xi,X2j «)} for the self-energy loop. Graph (b) vanishes since ipi and 
ipj cannot be on-shell simultaneously. 



Im{i 1 J 7< ^\xi,x 2 ;z)} =-Im{z 1 F > (x 1 , x 2 , x 3 ) - i 1 F < (x 1 , x 2 , z)} 



= Im 

2 



dH 



(2tt)< 



(+i gi )(-ig z )(-ig 3 )Dg(q + l)D^(q)D+(l)- 



(+igi)(-ig z )(-igs)D+(q + l)D^(q)D-(l)\s}, 



(27) 
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where (the result for) S coincides with eqn. (22) in the phenomenological scenario, while S — 1/2! includes 
an additional symmetrization factor in the toy model. This becomes 



Im{z 1 J r< ^\xi,x 2 ;z)} = 
1 f dH 



8{(q + I) 2 m$)5(l> - m^^-L^jl ± ft** + f^) S . (28) 



2 J (2tt) 2 ^ ' p > v a/ q 

This corresponds to the result for the self-energy contribution in the hierarchical limit in the top-down 
approach [5] if the equilibrium distribution functions are replaced with non-equilibrium ones. In the zero 
temperature limit this leads to the correct vacuum result. Thus, we have shown that (within the toy model) 
the CP-violating parameter e th obtained with help of thermal quantum field theory coincides with the one 
obtained in the top-down approach (in the approximately symmetric case) when one uses causal products 
instead of the conventional ones which assume type-1 external vertices. Furthermore, by comparing with 
the top-down result, we find that the thermal field theory result can be generalized to a (symmetric) non- 
equilibrium configuration for the toy model by the canonical replacement of the equilibrium distribution 
functions with the non-equilibrium ones: f eq —} f . 

For the phenomenological scenario we obtain in the limit of massless lepton and Higgs for the self-energy 
contribution (including a factor of 2, because the two components of the lepton doublet can propagate in 
the self-energy loop for a given transition) 

«,«, _ _ 1 ^Mj/Aft} M ' M > f dn 'n _..„.»..(, _ ,',;^ r i,w 



where E\p are again given by eqn. (15). In the zero temperature limit this reduces to the standard result 



S.vac = 1 v Im{(fet/ l )| j } MjMj 

8vr^ (hth)u Mf — Mf 



The complete CP-violating parameter is given by 



r th _ V.th S,th 



(30) 



where the vertex and self-energy contributions are given by eqns. (25) and (29) respectively. Therefore the 
overall dependence on the distribution functions (vertex and self-energy contribution) is given by 

i-f E : q +fk eq - 

In contrast to previous findings eqn. (16), this does not vanish in the limit when the Majorana neutrino 
decays at rest assuming massless £ and <j>. Therefore, it is qualitatively different from the conventional 
result. The new expression can lead to a significant enhancement of the CP-violating parameter, see fig. 8. 
Similar formulas can be derived for processes such as <j> —> N\l in the standard model (which can become 
relevant at higher temperatures) or for similar MSSM processes involving sneutrinos and sleptons. The size 
of the medium corrections depends primarily on the statistics of the particles in the loop, see fig. 9. 



V. BOLTZMANN EQUATIONS 



We can assume in addition that the structure of the Boltzmann equations for the phenomenological scenario 
is analogous to the one given in [4, 5] with appropriate quantum statistical factors for bosons and fermions 
respectively and appropriate symmetrization factors. This defines the full set of Boltzmann equations in- 
cluding medium corrections to the CP-violating parameter for the phenomenological scenario as derived 
above. With these modifications, the minimal network of quantum corrected Boltzmann equations for 
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<ei> M 
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FIG. 8: Temperature dependence of the CP-violating parameter in the Majorana neutrino decay relative to its vacuum 
value. Shown are the thermal average (e* ) / el ac (solid red line) and the values for various momentum modes e* /ej oc 
(dotted red lines) corresponding to |q| = T, —1 < sin(<5') < +1. For comparison we also show the conventional 
results ( e ^ h < conv ^ e va.c (-(j asne( j black line), where the leading effects cancel as described in the text. Equilibrium 
distribution functions for bosons and fermions with negligible chemical potentials are assumed. Note that the shown 
behavior can be modified if thermal masses are included, since the decay iVi — > i(j> (and the conjugate process) becomes 
kinematically forbidden if the thermal Higgs mass becomes too large. At even higher temperature the process <f> — > N\l 
becomes relevant instead [2]. 



thermal leptogenesis with hierarchical Majorana neutrino masses Mi <C M2, M3 takes the form (in ho- 
mogeneous and isotropic Friedman-Robertson-Walker space-time and not writing equations for the Higgs 
fields <fi, <fi which are considered to be in thermal equilibrium): 

L[/^](|k|) = C Nl ^[f N \f, /*](|k|) + C Nl ^[f N \f, /*](|k|) , (31a) 

L[f](\k\) = C^ Nl {f, f*, /^](|k|) , (31b) 

L{f%k\) = C^ Nl lfj*J Ni m\), (31c) 
where the Liouville operator is given by 

L[r\(x, k) = fc° (| - |k| H-^j f a (|k|) . (32) 

If the generated asymmetry is small, as we assume here, then f l f l and « In this case the 
CP-violating contributions to the right-hand side of eqn. (31a) cancel out and we obtain 

c Nl ^4f N \fj*m)+c Nl ^if N \f^f*m) 

x ({[i-/ik,]/ipi/r q i-/ik|[i--/ip|][i+/iti]} 

+ {[1 - /fi]/f P |/|t| - /ffl [1 - /f P |][l + /,*,]}) > 03) 
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Ci^ Nl [f,ftf Nl }(M =7i / duldU^(27r) 4 5(k+p- q )\M \ 2 Nl ^(k,p)[l-e 1 (\c l 



FIG. 9: Medium correction to the CP-violating parameters in the MSSM. The lines correspond to the thermal averages 
(e" 1 ) /el ac , and the shaded regions illustrate the momentum dependence of e\ h /ei ac for 0.5 < \q\/T < 4 and 5' = 0. 
Note also that in the weighted sum of Ni — > <j>£ and Ni — > <f>l processes the cancellation of the medium contributions, 
observed in the earlier publications, does not occur anymore. 

where, as usual, the tree-level amplitude for the Majorana neutrino decay is given by \A4o\ 2 Ni _± ei p (p, q) = 
l-Mo|jVi->.Z0 (P> l) = 2(/i^/i)n p ■ q. The collision terms for the (inverse) decay of the heavy particle into a 
(tcf) or a £<f> pair explicitly contain the CP-violating parameter e\ given in eqn. (30) but with the equilibrium 
distributions replaced by non-equilibrium ones f' eq — > f: 

Ci^ Nl [f, ft f Nl ](M) =~ J dlf p cfflf (2n) 4 5(k +p - q) \M \ 2 Ni ^ (k,p)[l + ei (|q|)] 

x{[i - + /,t|]/S - ffafU 1 - & - (34a) 

1 

2 „ 

X{[1 - /|k|][l + /|t,]/jjf - /|k|/|t|[l - /|q|]} ■ (34b) 

Note that the network of Boltzmann equations (31) should be understood in the generalized sense: the 
transition amplitudes differ from the usual perturbative matrix elements and do not have their symmetry 
properties as was noted in [4, 5]. The structure of the collision terms (34) differs from the conventional 
one. In particular, we did not include the processes £<f> 4-> £(j> explicitly, because the collision terms for 
the processes £cf> <-> Ni and £<f> «-» Ni do not suffer from the generation of an asymmetry in equilibrium. 
To obtain a consistent set of equations in the canonical bottom-up approach we would need to subtract the 
RIS part of the 5-matrix element for the processes £<fi o l(j>. Note, however, that it may be necessary to 
include the collision terms for £(f> O £<t> (derived in the top-down approach) in quantitative studies, because 
these can violate CP in general. Further scattering processes with top-quarks and gauge-bosons can also 
give relevant contributions. We note here that this result should be treated with care, because additional 
new effects could arise when the phenomenological scenario is investigated in the top-down approach. In 
addition, the applicability of the quasi-particle picture can not be tested in the framework of thermal field 
theory. In particular the results presented above will only apply in the hierarchical case [5]. The analysis of 
the resonant case requires the use of the Kadanoff-Baym formalism, which allows us to take into account 
the in-medium spectral properties of the mixing fields. 
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VI. CONCLUSIONS 



Inspired by a discrepancy between conventional results for the thermal corrections to the CP asymmetries 
in thermal leptogenesis and recent new results from non-equilibrium quantum field theory we have recon- 
sidered the calculation of the CP-violating parameters based on thermal quantum field theory. We find that, 
if causal products are used in the computation of the n-point functions, the results of both approaches can 
be brought into agreement in the framework of a toy model. We conclude that causal n-point functions 
must be used in the derivation of the CP-violating parameter in the phenomenological scenario as well. 
This leads to new expressions for the thermal corrections to the vertex and self-energy CP-violating param- 
eters. In contrast to the conventional results the thermal corrections do not vanish in the limit when the 
Majorana neutrino decays at rest assuming massless decay products. Therefore, it is qualitatively different 
from the conventional result and might give significant contributions to the generated baryon asymmetry. 
In the range from 0.1 to 10 of the dimensionless inverse temperature, thermal effects can enhance the CP- 
violating parameter by up to an order of magnitude. The asymmetry can be computed using the minimal set 
of Boltzmann equations for leptogenesis in SM+3vr presented here, which are analogous to the equations 
which have been derived earlier in the framework of the toy model. These take into account decays and 
inverse decays and include all quantum statistical factors in a way which guarantees that no asymmetry is 
generated in equilibrium. They can be applied in the case of non-degenerate Majorana neutrino masses. For 
a detailed phenomenological analysis it will be necessary to take into account further thermal effects such 
as thermal masses and resummed thermal propagators as well as additional CP-violating processes which 
exist in phenomenological scenarios. 
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